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We introduce a two dimensional nonlinear XY model with a second order phase transition driven 
by spin waves, together with a first order phase transition in the bond variables between two "bond 
ordered phases" , one with local ferromagnetic order and another with local anti-ferromagnetic order. 
We also prove that at the transition temperature the bond-ordered phases coexist with a disordered 
phase as predicted by Domany, Schick and Swendsen-. This last result generalizes the result of 
Shlosman and van Enter^. We argue that these phenomena are quite general and should occur for 
a large class of potentials. 

PACS numbers: 64.60.Cn,75.10.Hk 



I. INTRODUCTION AND MAIN RESULT 

A usual feature in the Statistical Mechanics is the coex- 
istence of several phases at low temperature. Typically, 
the number of coexisting pure phases decreases with in- 
creasing temperature up to some critical temperature f3c 
at which a single disordered phase appears. This g-state 
Potts model for sufhciently large q behaves slightly differ- 
ently since it has q ordered phases below the transition 
temperature, but at the transition it has q + 1 phases, 
including a disordered phaseS^i. Thus, the Potts model 
demonstrates that order and disorder can coexist. 

There is a continued interest in exploring the variety 
of ordering phenomena occurring in models of classical 
continuous XY spins^i^. In the present work, one of our 
motivations was to look for the coexistence of ordered 
and disordered phases in such models (as opposed to the 
previously studied Potts model, which is a discrete spin 
model) . 

An interesting situation of coexisting order and dis- 
order was found to occur in a strongly non-linear XY 
modeli*^. For that model it was recently proved rigor- 
ously by Shlosman and van Enter in Ref. |2 that there is 
some temperature at which a first order phase transition 
occurs between a (ferromagnetic) bond ordered phase, 
which means that the nearest neighbors are close, and 
a bond disordered phase, which means that the nearest 
neighbors are uncorrelated. This behavior resembles the 
Potts model with a magnetic field, which shows a unique 
ordered phase at low temperature, and two coexisting 
phases at the transition. The main difference is that the 
XY models has a continuous symmetry (0(2)), as op- 
posed to Potts model where the symmetry is discrete. 
This continuous symmetry cannot be spontaneously bro- 
ken by the Mermin- Wagner Theorem, which implies any 
long-range order is not accompanied by a non- vanishing 
magnetization. 

In this paper we introduce a new family of 0(2)-models 



defined on the two-dimensional square lattice, some of 
which exhibit a second order phase transition at suffi- 
ciently low temperature that coexists with a first order 
phase transition. Our family of models depends on five 
parameters: three positive integers, p, m, and n, and two 
non-negative coupling constants, J, and K: 



H = 2J^ [cosy(v5i 



2p 



K 



/ , cos -((^j - lpk) 
[[i,fe]] 

(1) 

where ipi g [0, 27r), and [., .] denotes nearest neighbors 
and [[.,.]] denotes diagonal neighbors. The model con- 
sidered in Refs Q and 2 is the special case K — Q and 
TO = 1. As we will only consider X > 0, we can fix K by 
rescaling the temperature; we pick K — 1/4. 

We prove that, for suitable values of the parameters, 
these models have both a first and a second order phase 
transition. Our main result is the following theorem, 
in which we demonstrate the occurrence of both phase 
transitions by means of estimates on on nearest neighbor 
and long range correlations. The correlation functions 
with specified boundary conditions (b.c.) are denoted 
by (•)'' '^'. We consider three types of boundary condi- 
tions: 1) /, for ferromagnetic, indicating parallel nearest 
neighbor spins, 2) a, for antiferromagnetic, which favors 
antiparallel nearest neighbor spins, and 3) d, for disor- 
dered, denoting the third phase, which appears at the 
transition temperature. The precise meaning of the var- 
ious boundary conditions is clarified later. 

Theorem. For sufficiently large p, m — 1, and n = 

2, the models with Hamiltonian 0) have the following 
properties: 

I) There exists an inverse temperature (3t, and func- 
tions €°''^'{(3,p) < 1, and i]^''^'{(3,p) < 1, such that: 

a) for P > Pt, there is a first ordered phase transition 
between the bond ferromagnetic phase and the bond anti- 



ferromagnetic phase: 

{cosiip, - ipj))f > l-ef{f3,p) (2) 
(cos(^, - ^,))'^ < -l + e«(/3,p) (3) 

for nearest neighbor pairs [i,j]. 

b) at Pt there is a first order phase transition between 
the bond ordered ferromagnetic phase, the bond anti- 
ferromagnetic phase and the disordered phase: 

(cos((p, -(^j)>^ > ^-vHp) (4) 
(cos((^,-^,))" < -l + v^'ip) (5) 
|(cos((p.-^,))'^| < v'ip) (6) 

II) At sufficiently low temperature there is a second order 
phase transition to phase with power law decay of the two 
point correlation: 

C 1 

— < {cos{ipo - (Pn)) < — 1 (7) 

where (3' > and /3'(/3) oo when /3 — > cxi, and N 
denotes any site at distance N from the origin. 

We interprete our result as follows: two phases coexist 
at low temperature, one with a local ferromagnetic or- 
der, and one with a local anti-ferromagnetic order, but in 
large distance spin waves destroy both the long range fer- 
romagnetic order and the long range anti-ferromagnetic 
order leading to a second order phase transition at some 
lower temperature. The situation has to be opposed to 
the rotator model in two dimensions, which exhibits a 
second order phase transition at low temperature, but no 
first order phase transitior&iS. Notice that the Lebowitz 
inequalitiesii which were used in Ref. to prove the 
uniqueness of the correlation functions in the two dimen- 
sional rotator model are not valid in our model. 

In the next section we describe the different phases 
of our model more precisely and introduce associated re- 
stricted ensembles used for the proof of the theorem. The 
proof itself is given in a separate section. 

II. THE BOND-ORDERED PHASES. 
RESTRICTED ENSEMBLES 

We want to prove that below a temperature two 
two order phases coexist: one in which ferromagnetically 
ordered bonds dominate and another in which antiferro- 
magnetically ordered bonds dominate. We also want to 
show that at the first-order transition temperature, /3t, 
three phases coexist: the two ordered phases plus a dis- 
ordered phase. We will give the complete proof of this 
last statement, from which it will be then be clear how 
to prove the first. First, we give some definitions and 
describe the ground states of the model. 

For the purpose of describing the ground states, a 
plaquette is called ferromagnetic if ipi = (pj for each 
nearest neighbor bond, [i,j], of the plaquette, and anti- 
ferromagnetic if ipi — ifj = TT for each bond. We will 
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FIG. 1: Three examples of plaquettes. Thick lines denote fer- 
romagnetic bonds (nearly parallel spins), thin lines denote an- 
tiferromagnetic bonds (nearly antiparallel spins) and a zigzag 
line is used to denote all other cases, which we call disordered 
bonds. The first two plaquettes are the ferro and antiferro- 
magnetic ground state plaquettes and will be labeled //// 
and aaaa, respectively. The third plaquette is neither ferro- 
magnetic, antiferromagnetic, nor disordered and is labeled by 
the sequence adaf. 



extend these definitions to finite temperature a couple of 
paragraphs further down. 

The Hamiltonian ijlp-^'^} has two families of ground 
states, each parametrized by an angle: 

• the ferromagnetic ground states, in which all plaque- 
ttes are ferromagnetic; 

• the antiferromagnetic ground states, which have only 
antiferromagnetic plaquettes. 

It is not hard to see that these are the only ground 
states and it is easy to specify boundary conditions that 
select one of them. 

Next, we define families of configurations that are close 
to one the ground state configurations and that will carry 
most of the weight of the equilibrium states for an appro- 
priate range of temperatures. Restricted ensembles and 
the corresponding partition functions can then be de- 
fined by summing over all configurations in each of these 
families. For each bond b = [i,j] three characteristic 
functions: 

Here again, and in the sequel, / stands for ferromagentic, 
a for antiferromagnetic, and d for disordered. 

The specification of the bond is ferromagnetic if X'^' [•] — 
1, anti-ferromagnetic if x'^H-] ~ ^^'^ disordered if 
X^[-] = 1- In the figures we will indicate a ferromag- 
netic bond by a thick line, an antiferromagnetic bond by 
a thin line, and a disordered bond by a zigzag line. The 
specification of the four bonds of a plaquette will be de- 
noted by stuv, where s, t,u,v S {/, a, d}, and where the 
bonds are ordered in clockwise fashion starting from the 
top horizontal bond. As before, a ferromagnetic plaque- 
tte is one with specification ////, an antiferromagnetic 
plaquette has aaaa, etc. (See Figure QJ. 
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configuration there corresponds an antiferromagnetic one 
obtained by adding tt to each spin on the sites of A. 

The restricted partition functions Sv(/), Sy(a), and 
E!v((i), are the partition functions with the configuration 
sum restricted to Ry, Ry, and Ry, respectively. 

The one-to-one correspondence between Ry and Ry, 
imphes that Ey = Ey. 

A contour T is the set of configurations, which are 
equivalent to a precontour's configuration 7. A contour 
r is specified by: 1) a geometric support S'(r), which is 
a connected set of plaquettes; 2) a family of equivalent 
configurations defined on S(r). 



FIG. 2: The universal contour U{qfaf) obtained from the pla- 
quette specification [Zl by repeated reflections. 



III. THE PROOF 



FIG. 3: The universal contour U{adaa) obtained from the 
plaquette specification HI by repeated reflections. 



For a general configuration, we say that 7 is a precon- 
tour if it is is a maximal connected set of plaquettes that 
are neither ferromagnetic nor antiferromagnetic nor dis- 
ordered. A precontour is composed of different plaquettes 
i.e. plaquettes containing two or three different bonds: 
either ferromagnetic, or anti-ferromagnetic or disordered. 
A precontour can be thick, and different precontours can 
have the same support. Notice that below the pre- 
contours contain typically disordered plaquettes. 

The configurations of two plaquettes are called equiv- 
alent if the specification of each plaquette is the same. 
An equivalent class of configurations in is a restricted 
ensemble. 

We define three kinds of restricted ensembles: 1) the 
restricted ferromagnetic ensemble Ry built with equiva- 
lent ferromagnetic plaquettes in V; 2) the restricted an- 
tiferromagnetic ensemble Ry built with equivalent anti- 
ferromagnetic plaquettes; 3) the restricted disordered en- 
semble Ry built with equivalent disordered plaquettes. 

There is a one to one correspondence between the re- 
stricted ensembles Ry and Ry. To see this, consider the 
decomposition of the lattice into the two standard sub- 
lattices A and B, where the nearest neighbors of each 
site in B (A) are in A (B). The, with each ferromagnetic 



Proof of I. a and I.b. The existence of a first order 
phase transition^ is a consequence of the three conditions: 

a) both the probability of an ordered plaquette of an 
antiferromagnetic plaquette are small at high tempera- 
ture, 

b) the probability of a disordered plaquette is small 
at low temperature, c) the probability of the contours 
are small at every temperature. The proof of La) and 
I.b) will be achieved via chess-board estimates, which 
require that the Hamiltonian is reflection positiv^^ti^'^. 
We expand the Hamiltonian. 

H = J A„i cos2m{(pi — ipj) 

[i,j],m=l,2,...q) 



+ J ^^^ifi - Vk) 



(11) 



[[i,fc]]ev 



Here, the Am are positive integers. Obviously the 
Hamiltonian is reflection positive with respect to 
the axis of the lattice. For simplicity we choose J = 1. 

We prove directly part I.b), because I. a) follows by 
symmetry. Peierls bounds via chessboard estimates 
are obtained using universal contours, which are ob- 
tained from a plaquette configuration stuv, by reflections 
through the axes of the lattice. See Figures 2 and 3 for ex- 
amples of universal configurations U{stuv), which should 
be thought of as equivalence classes of configurations. 

The probability of a plaquette configuration stuv is 
given via the following formula, which follows from suc- 
cessive applications of the Schwartz inequality iS^ii. 



(stuv) >v< 



Ev[U{stuv)] 



(7) 



First we need to prove a lower bound for :^v, which is 
similar to the one used in Ref. 



5v > 



\v\ 



e^pP\V\[^-0{^)] + 



AC 



We briefly recall the argument leading to H12|) . 



(12) 



FIG. 4: The universal contour U{adad) obtained from the 
plaquette specification s_j by repeated reflections. 



FIG. 6: The universal contour U{dadd) obtained from the 
plaquette specification 'C] by repeated reflections. 




get the upper bound for Sy[/7(a/a/)] by noticing that the 
energies of the diagonal potentials are zero: 



EvPiafaf)] < 



a) We get for P > Po- 



\v\ 



exp[/3|y| 



2 + 0( 



C2- 




FIG. 5: The universal contour U{aadd) obtained from the 
plaquette specification O by repeated reflections. 



The first term is obtained by integrating over all the 



configurations such that \tpi 



< 



for all i £ V, the 



second one by integrating over all the configurations over 



> 



for all 



one Neel sublattice which satisfy \tpi 
nearest neighbors. 

Secondly we derive an upper bound for the quantities 
Ev[U{stuv)] for the different restricted ensembles P{D). 
The proofs of the conditions a) and b) are straightfor- 
ward, so we focus on the proof of c), which is to prove 
that the probabilities of the various plaquettes occurring 
in contours are small at every temperature. Let [3q be 
defined by: 



(13) 



We have to estimate the probability of the contours 
with different specifications. 

• Two opposite bonds with = 1 and two opposite 
bonds with x° = 1, i.e., a plaquette with specification 
afaf. The universal contour U{afaf) is built from ferro- 
magnetic vertical and horizontal bonds with even coordi- 
nates, and from anti-ferromagnetic vertical and horizon- 
tal bonds with odd coordinates as shown in Figure[21 We 



(n) 



< 



\v\ 



cvp. 



xexp/?|y|(|-0(^)) + 



1- 4C 



< exp-/3^i-0(ij) 



< 



1 



AC 
VP. 



(14) 



b) For /3 < /3o we obtain: 



(n> < 



pc\i-f,) 



[§+0(^)1 



cVp 



1 _ 4C 

VP. 



< 



cVp 



Vp. 



In both cases the expectations can be made small for p 
large enough compared to C. It is obvious that the upper 
bounds for the expressions of (IZI) computed at low and 
at high temperature have the same dependence on C and 
p. This holds in general and therefore we will compute 
only one case for the other specifications. 

• One disordered bond and three ordered bonds. The 
nature of the ordered bonds, / or a, produces only minor 
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changes. Therefore, let us focus on one instances, say 
the universal contour U{adaa) which is constructed from 
ordered horizontal bonds and even vertical bonds, and 
from disordered vertical odd bonds, as shown in Figure 
O For the upper bound for 'E,v[U{adaa)], we notice that 
the energy of the odd vertical bonds is zero: 



^v[U{adaa)] < 



1 



|V| 



which leads to the upper bound 
1 



expp\V\[2 + 0{-^)], 



(□) < 



cVp. 



(15) 



• Two opposite disordered bond and two opposite or- 
dered bonds. The universal contour U (adad) is built from 
ordered horizontal bonds, and from disordered vertical 
bonds as shown in Figure ^ We get the upper bound 
for the 'B.v[U{adad)] by noticing that the energy of the 
vertical bonds is zero: 



Ev[U{adad)] < 



1 



We get the upper bound: 

(n)< ^ 



expP\V\{^+Oi-^) (16) 



i(f-o[^l) 



1 + 



4C 
VP 



(17) 



• Two adjacent disordered bonds and two adjacent or- 
dered bonds. The universal contour U (aadd) is built 
from ordered horizontal even bonds and even vertical 
bonds, and from disordered vertical odd and the hori- 
zontal bonds as shown in Figure We get the upper 
bound for ^v[U{aadd)] the restricted partition function 
by noticing that the energy of the odd vertical and hori- 
zontal bonds are zero, and that there are entropy contri- 
butions from the vertices with odd coordinates: 



Ev[U{aadd)] < 



1 



\\v\ 



cVp\ 

We get the upper bound: 



exp[/3|y| 



^0(i,) 



1 



CVp 



-o[^]) 



v^. 



(18) 



(19) 



• Three disordered bond and one ordered bonds. The 
universal contour U (dadd) is built from disordered hori- 
zontal bonds and even vertical bonds, and from ordered 
vertical odd bonds as shown in Figure El We get the 
upper bound for ^v[U{dadd)] by noticing that the en- 
ergy of the horizontal and the odd vertical and bonds is 
zero, and that there are entropy contributions from the 
vertices with even ordinates: 



Ev[U{dadd)] < 



1 



CVp 



exp/?|F|[l + 0(i^)] (20) 



We get the upper bound: 



< 



CVp. 



4C' 
VP. 



(21) 



For p ^ C ^ 1, the probability of each specification of 
a contour's plaquette can be made small enough. Finally 
we prove the Peierls condition: first we sum over all the 
specifications of the contours with a fixed support using 
the estimates (^51, 03, OHl, and 121]); next we 

use the Koenigsberg's lemma to get an upper bound on 
the entropy for the contours. This concludes our proof 
of part I of the Theorem. 

Proof of II) The upper bound on the two points cor- 
relation function is proved via a spin wave argument, 
the lower bound is deduced from the deep result of 
Frohlich and Spencer for the two dimensional classical 
XY model^^. We now provide the details. 

II. a) The lower bound 

We use the Ginibre inequalitiesi^, which are valid for 
ferromagnetic Hamiltonians. 

{cos{mi[ipi - ipj]) X cos{m2[(pk - fi])}"^ 
> (cos(mi[<y9, - (pj]))"^ {cos{m2[(pk - V'/]))^^' (22) 

We first start by removing all the terms contained in 
the Hamiltonian with the exception of the diagonal in- 
teractions, then we are left with two independent rotator 
models, one on each of the sublattices. 

First, suppose i and j are on the same sublattice. 
Then we use the remaining interactions on that sub- 
lattice. We are left with the Hamiltonian: Hy =: 
E{[[i,fc]]ev} cos([(^fe - (pi]), and we have 



(cos((^o - Vn))""-' > (cos(v3o - Vn))""^ 

Next we apply the result of Frohlich and Spencer^'^ , The- 
orem C, for the two dimensional rotator model to get a 
lower bound. 



(cos((^o - Vn))"'^ > 



C 



(23) 



where fi' > and /3'(/3) — > oo when /3 — > oo. 

If i and j belong to different sublattices, the same lower 
bound of 123|l holds with a slightly smaller constant C, 
again by using the Ginibre inequality (|22|) with j = I and 
k and j replaced by a nearest neighbor of j . 

II. b) The upper bound. The upper bound for the 
two points correlation function is obtained by using a well 
known spin wave argument. We use a nice argument that 
we learned from H. Kunaii. It starts by defining a se- 
quence of squares (or hypercubes in higher dimensions) 
and label them Bi, B2, . . ., as shown in Figure [7| Each 
Bi is centered at the origin and its boundaries can be 
taken to intersect nearest neighbor bonds. We define the 
shells Ci = -Bj+i \ Bi, each shell C; contains a square 
Ri of the lattice, the nearest vertices contained in Ri are 
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FIG. 7: The nested square boxes Bi used in the argument 
for the upper bound II. b. are shown in dashed lines. The 
solid bonds are the new strongly ferromagnetic interactions 
introduced in the proof. 



labelled by a couple [I, where denotes the coordi- 
nates along Ri . We want to identify the continuous spins 
contained in each shell Ri by using Ginibrc's inequality 
(|22|l . To do so, we add to the Hamiltonian Hy the ferro- 
magnetic interactions: Y.\Zi Y.i{i) Ji cos(<^i+i(i) - ipi(i)), 
where the J; is are positive. The correlations are increas- 
ing by Ginibre's inequality. Next we let Ji — oo, the 
correlations are again increasing: 



(cos(i^o - V'n))v 



< 



1 



dipi--- / TT X 

exp [4jJ AmCos2m{ipj — (fj^i) 

m— l,2,...g 



Where Z'f^ is the normalization factor. The LHS is 
the two points correlation function corresponding to a 
one dimensional model, each term of the product can be 
computed, we get the following upper bound: 



IV. CONCLUSION 



The coexistence between second order and first order 
phase transitions at the same temperature should be a 
fairly general phenomena in two-dimensional models with 
a continuous symmetry. For example it occurs in several 
non linear two dimensional XY models defined by (1) for 
other values of the pair m, n. In a more general setting, 
we have to consider potentials with continuous symmetry, 
which are peaked and with large flat parts to produce en- 
tropy. The problem will be to characterize the criticality, 
which in our case uses^'^, a generalization of the work of 
Aizenman could be an alternativei^. We expect that the 
same situation should hold for some quantum models. 
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